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ABSTRACT 

The cosmic infrared background (CIB) is a sensitive measure of the star formation history. But 
this background is overwhelmed by foregrounds, which bias the CIB mean flux and auto correlation 
measurement severely. Since dominant foregrounds do not correlate with the large scale structure, the 
cross correlation of CIB with galaxies is free of such foregrounds and presents as an unbiased measure 
of the star formation history. In such cross correlation measurement, one can utilize all frequency 
information, instead of being limited to several narrow frequency windows. This allows the measurement 
of CIB based on integrated intensity, whose theoretical prediction is based on energy conservation, thus is 
fairly model independent and robust. The redshift information of CIB sources can be directly recovered 
with the aid of galaxy photometric rcdshifts. In the far IR region, correlation signal is around 10% at 
degree scale. Combining FIRAS and SDSS, the cross correlation can be measured with 20% accuracy 
(statistical and systematics errors) at I ~ 40. Such measurement would allow the measurement of star 
formation rate with comparable accuracies to z ~ 1. Future CIB experiments with 1° resolution are 
able to measure the star formation rate with 10% accuracy to z ~ 1.5. Secondary CMB anisotropies 
arising from the large scale structure correlate with CIB. The correlation strength is about 0.5/iK at 
~ 10° (the integrated Sachs- Wolfe effect) and ~ — 0.3/xK at ~ 1° (the Sunyaev Zeldovich effect, Rayleigh- 
Jeans regime). FIRAS+WMAP would allow the measurement of the ISW effect with > 5er confidence. 
Its detection helps to constrain the amount of dark energy and far IR emissivity bias. A future CIB 
experiment with better angular resolution would allow the measurement of CIB-SZ cross correlation and 
help to understand the effect of SN feedback on the SZ effect. 

Subject headings: cosmology: large scale structure; infrared: theory-diffuse background; star: formation 
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1. INTRODUCTION 

Detections of the cosmic infrared background (CIB) (re- 
fer to Hauser & Dwek (2001) for a recent review) open a 
new window for the study of star formation history. The 
main CIB component associates with star formation. Dust 
in extragalactic galaxies absorbs UV photons and reemits 
in far infrared region. Cosmic expansion shifts light emit- 
ted by low mass population I/II stars or massive popu- 
lation III stars into near infrared region. We denote this 
CIB originated from star formation as SCIB, which pro- 
vides a fairly unbiased and statistically robust measure of 
the star formation history (SFH). But the precision CIB 
measurement is challenged by overwhelming foregrounds. 
Interplanetary dust emission dominates over the CIB by 
an order of magnitude at A < lOO^iin. Galactic interstel- 
lar dust emission peaks at A ~ 100/im with an intensity at 
least comparable to the CIB. Bright galactic sources such 
as stars and faint galactic sources have intensities at least 
several times larger than the CIB in near IR region. The 
CMB dominates over the A > 400/im far IR region. 

To reliably separate the CIB and these dominant fore- 
grounds, robust understanding of foreground spacial and 
spectral distribution is required. The removal of these 
foregrounds has enabled successful CIB detections. But 
the residual foregrounds still cause large dispersions in the 
current CIB intensity measurements (see table 1 of Hauser 
& Dwek (2001) for a review up to 2001 and Wright (2003) 
and reference therein for a recent compilation of COBE 
data). 

One way to bypass foregrounds is to investigate the 



CIB fluctuations, as proposed by Bond, Carr, & Hogan 
(1986); Haiman & Knox (2000); Knox, Cooray, Eisenstein, 
& Haiman (2001). This fluctuation approach has several 
advantages. (1) The CIB roughly follows star forming 
galaxy (SFG) distribution, so one expects ~ 10% fluctua- 
tions at degree scale, which are about 4 orders of magni- 
tude larger than the CMB fluctuations. (2) Fluctuations 
of galactic foregrounds concentrate on large angular scales, 
for example, the power spectrum of the galactic interstellar 
medium emission scales as Q cx l~ 3 (Wright 1998). Thus 
at sufficiently small angular scales, The CIB fluctuations 
exceed these foreground fluctuations. 

But correlations of CIB foregrounds bias this approach. 
Fluctuations of Galactic interstellar dust emission dom- 
inate over the CIB intrinsic fluctuations in a large sky 
fraction and large frequency and angular ranges (Knox, 
Cooray, Eisenstein, & Haiman 2001). The CIB auto corre- 
lation measurement in such regions requires an accurate re- 
moval of these foreground correlations. The residual fore- 
ground correlations bias the intrinsic CIB auto correlation 
measurement and the following extraction of SFH. This 
problem is most severe at large angular scale (I < 200) 
or short wavelength (A < 300/xm) and may affect the es- 
timation of Knox, Cooray, Eisenstein, & Haiman (2001) 
significantly. To minimize such bias, the auto correlation 
measurement has to be limited to finite sky regions, fre- 
quency bands and angular scales. These limitations in- 
crease the sample variance. 

For the purpose of extracting star formation rate (SFR), 
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various non-stellar CIB, such as AGN CIB 1 , introduce ex- 
tra systematic errors. AGNs have similar clustering prop- 
erty as galaxies and similar thermal emission features as 
extragalactic dust, so it is hard to separate AGN CIB 
and SCIB by the mean flux and auto correlation measure- 
ments, even utilizing multi-frequency information. Since 
AGN may contribute ~ 10% to the CIB intensities (Hauser 
& Dwek (2001) and reference therein), it may introduce 
<~ 10% systematic error to the SFR measurement. 

The ultimate challenge of the auto correlation approach 
is to recover the redshift information of CIB sources. 
This can be done by multi-band CIB correlations (Knox, 
Cooray, Eisenstein, & Haiman 2001). But such approach 
requires detailed understanding of dust composition and 
size distribution and UV sources distribution, such mod- 
eling is difficult and uncertainty it introduces is hard to 
quantify. 

With all these concerns, new CIB analysis method is 
demanded. Dominant CIB foregrounds such as galactic 
foregrounds and primary CMB do not correlate with extra- 
galactic galaxies. This suggests that crossing correlating 
the CIB with galaxies may alleviate foreground contami- 
nations. Residual correlations originated from other CIB 
and extragalactic contaminations arc negligible compar- 
ing to the SCIB-galaxy correlation. Though AGN CIB- 
galaxy cross correlation still exists in such approach and 
could bias the result, its effect is at most ~ 10%. Fur- 
thermore, one can utilize the difference of AGN redshift 
distribution versus that of SFGs to at least partially sep- 
arate the AGN contribution. Thus, the cross correlation 
measurement of the CIB and galaxies is fairly unbiased. 
With the aid of galaxy photometric redshift information, 
one does not need to rely on frequency information to re- 
cover redshift information of CIB sources. We only need 
to model the integrated CIB intensity, whose prediction is 
directly based on energy conservation and thus fairly ro- 
bust. So the cross correlation of CIB with galaxies would 
avoid most problems in the mean flux and auto correlation 
approaches. 

Since dust emission in far IR region is directly related to 
SFR, throughout this paper, we focus on the far IR region. 
Similar method can be applied to near IR region, which 
will constrain SFR at low mass end. In §2, we estimate the 
auto and cross correlation power spectra of SCIB based on 
the integrated CIB intensity. Our model can be applied 
to CIB-lensing cross correlation, as Song, Cooray, Knox, 
& Zaldarriaga (2003) did. But weak lensing surveys are 
limited by finite sky coverage while redshift information 
of lenses is entangled, so we postpone such calculation. 
We investigate the observational feasibility of such cross 
correlation by estimating the systematic errors (§3.1) and 
the statistical errors (§3.2). Since the SCIB traces the 
large scale structure, it correlates with secondary CMB 
anisotropies such as the integrated Sachs- Wolfe (ISW) ef- 
fect and the thermal Sunyaev-Zeldovich (SZ) effect. The 
cross correlation of the CIB and the CMB provides an- 
other way to constrain dark energy and the large scale 
structure. We predict the strength of such correlation in 
§4. Our predictions are based on the WMAP-alone cos- 



mology with f2 m = 0.268, Qa = 1 — Q m , 



0.044, 



cr 8 = 0.84 and h = 0.71 (Spergel et al. 2003) 2 . 

2. INTRINSIC COSMIC INFRARED BACKGROUND 
FLUCTUATIONS 

We model the SCIB based on energy conservation. Ob- 
servations suggest that /absption ~ 90% UV light is ab- 
sorbed by dust (Massarotti, Iovino, & Buzzoni (2001) and 
reference therein) and reemitted in far IR. Since UV emis- 
sion is dominated by short- life O/B stars, UV emissivity 
puv is directly related to SFR by 

Coefficient a„ is determined by the initial mass func- 
tion (IMF). It strongly depends on the shape of IMF 
at high mass end. At 1500 A, for the Salpeter IMF, 
a v ~ 8.0 x 10~ 25 nW Hz _1 (Madau, Pozzetti, & Dickinson 
1998). For the Scalo IMF which has much less massive 
stars, a v is about twice smaller. Since the Scalo IMF pro- 
duces too red integrated galaxy spectra (Lilly, Le Fevre, 
Hammer, & Crampton 1996) and too low CIB intensity, 
we consider only the Salpeter IMF in this paper. For the 
Salpeter IMF, a v has only a weak dependence on v. For 
example, at 2800 A, a v ~ 7.9 x 10~ 25 nW Hz" 1 . Then 
the integrated dust FIR emissivity is given by 

-■ - SFR 
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SCIB intensity at direction n is 

T(fl s _ [ jdQ- + Sj{x,n)) 
() -J 4tt(1 + z) 2 
Here, Sj is the fractional fluctuation of the integrated 
IR emissivity. \ lS ^ ne comoving distance. SFR at 
z < 1.2 is relatively robustly measured and we adopt it as 
SFR(z) = lO- 21 exp(t/2.6Gyr)M yr- 1 Mpc" 3 (Hippelein 
et al. 2003). Here, t is the look-back time. For z > 1.2, we 
consider three SFR: 

t-t(z = 1.2) 



SFR(z = 1.2)exp(- 
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For the Salpeter IMF, we obtain I ~ 16, 20, 30 nW m~ 2 sr" 1 , 
respectively. These results are consistent with observa- 
tions, but due to the large dispersion in the observational 
results (see Hauser & Dwek (2001); Wright (2003) for 
reviews), no specific conclusion can be drawn. 

We calculate the CIB fluctuations by the Limber's equa- 
tion: 

2 ^(Lz)d X , (5) 
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Here, A 2 and A 2 G are the variances (power spec- 
tra) of Sj and its cross correlation with galaxy num- 
ber overdensity Sq, respectively. dn/dz(z) is the 



1 We always refer AGN CIB as the CIB contributed by AGN accretion. 

2 For simplicity, we adopt a slightly different Ha from the best fit value 1.02 — Cl m . This change has negligible effect on our results. 
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galaxy number distribution function in a given galaxy 
survey. The galaxy surface density is defined as 
Eg = J dn/dz(l + 5o)dz. We adopt the parametric 
form dn/dz = 3z 2 /2/(z m /f .412) 3 exp(-(1.412z/z m ) 1 - 5 ) 
(Baugh & Efstathiou 1993). z m is the median redshift 
of the galaxy distribution. For SDSS, We z m ~ 0.5 
(Dodclson ct al. 2002). The above equations hold at 
angular scales < 30° or I > 10. For smaller /, one 
has to calculate the spherical Bessel integral. Since the 
small I modes are dominated by cosmic variance and/or 
foregrounds, whose power concentrates on large angu- 
lar scales, we neglect the calculation for such small I. 

10 100 1000 10 4 




Fig. 1. — The CIB power spectra. The top panel plots the CIB 
auto correlation power spectra while the bottom panel plots the 
cross correlation power spectra with SDSS galaxies. We assume 
three SFR: SFR1 (dash lines), SFR2 (dot lines) and SFR3 (solid 
lines) . The change of amplitude reflects the change of redshift where 
the dominant correlation signal comes from. 

A 2 and A 2 G are generally unknown. Since dust associates 
with SFGs, averaging over many galaxies, jd should trace 
SFG distribution, which in turn traces dark matter distri- 
bution. A 10 2 sky area contains more than 100 SFGs, so 
one expects that the constant bias assumption 5j = bj Sum 
applies at > 3 scale. In principle, bj is an observable and 
can be inferred from the combined CIB mean flux and cor- 
relation measurements (§3.2) or from the CIB-ISW cross 
correlation measurement (§4). For simplicity, we assume 
bj = &sfg = 3 (Haiman & Knox 2000). This SFG bias 
is inferred from Lyman break galaxies (Giavalisco et al. 
1998). The dark matter power spectrum is calculated by 
the code of Smith et al. (2003). 

The fractional fluctuations at degree scales are about 
10% at degree scale, which translate to intensity fluctua- 
tions ~ 2 nW m~ 2 sr _1 . These fluctuations are consistent 
with both theoretical predictions (Haiman & Knox 2000; 
Knox, Cooray, Eisenstein, & Haiman 2001) and observa- 
tions (Hauser & Dwek 2001; Miville-Deschenes, Lagache, 
& Puget 2002). 



3. ERROR ESTIMATION 

The above prediction assumes that the observational 
wavelength band contains all extragalactic dust emission. 
The requirement on the wavelength band can be esti- 
mated as follows. Extragalactic dust emission can be de- 
scribed by a gray-body spectrum I v ex v$ B V (T) with typ- 
ical T ~ 20 K (see, e.g., Haiman & Knox (2000)) and 
P = 2 (Draine & Lee 1984). Here, B V (T) is the Planck 
function. vl v peaks at ~ 2.5T/20K THz (or, equivalently, 
20K/T120/im). More than 96% emission comes from the 
range [0.8,5]T/20K ([60, 375]20K/7>m). Since star for- 
mation rate was more than an order of magnitude higher 
at z ~ 1 than at present, considering the redshift effect, an 
integral range [0.3,3] THz ([100, 1000]^m) should contain 
almost all contribution from extragalactic dust emission. 
Indeed, from the fitted FIRAS far infrared background 
(FIRB) spectrum (Fixsen et al. 1998) 

I„ = 1.3±0.4xl0- 5 (— ^-)°- 64±0 - 12 B,(18.5±1.2K), (7) 
3 1 Hz 

this fraction is > 92%. This frequncy range may under- 
estimate low redshift contribution by 10%, but it contains 
> 95% extragalactic dust emission from z > 0.6 and is 
sufficient for our extraction of SFR at z ~ 1. 

FIRAS satisfies this wavelength band requirement, so 
we take it as our target to estimate both the statistical 
and systematic errors of the cross correlation measure- 
ment. The assumption that [0.3,3] THz contains all dust 
emission would introduce several percent errors at low red- 
shift, which are minor comparing to other systematic er- 
rors (§3.1). 

3.1. Systematics errors 

Non-stellar CIB sources, if correlate with large scale 
structure, bias SFR measured from the CIB-galaxy cross 
correlation. In this section, we estimate the systematic 
errors caused by such sources. We find that the dominant 
systematic error comes from AGN CIB, which may con- 
tribute 10% to the cross correlation. We show that AGN 
CIB can be at least partially subtracted. In the worst case 
that these sources are not subtractive, they can introduce 
at most 10% systematic error to the measured jd, which 
is still quite accurate comparing to current SFR measure- 
ment (see, e.g. Hippelein et al. (2003) for recent data 
compilation) . 

3.1.1. AGN 

AGN may contribute ~ 10% to the CIB flux through 
the thermal emission of dust heated by the central black 
hole (Hauser & Dwek (2001) and reference therein). This 
non-stellar AGN contribution would affect the determina- 
tion of SFR by a comparable fraction. AGNs have similar 
clustering property as galaxies and similar thermal emis- 
sion feature as extragalactic dust, so it is hard to separate 
AGN CIB and SCIB by the mean flux and auto correla- 
tion measurements, even with the aid of multi-frequency 
information. Many high redshift AGNs present as lumi- 
nous infrared galaxies and dominate Lm > 10 11 L Q object 
catalog (Sanders & Mirabel 1996). So point source iden- 
tification, at least at high luminosity end, may be feasible 
and help to understand the contribution of AGNs. But the 
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distinguishing of energy sources of such IR emission, grav- 
itational potential energy or nuclear energy, is not very 
clear ((Hauser & Dwck 2001) and reference therein). 

AGN contamination can be partially cleaned by cross 
correlation with galaxies. 

• Dust heated by starbursts has temperature T ~ 
20K while dust associated with AGN IR emission 
has temperature T ~ 60-100 K (Haas et al. 1998). 
Since the observed Far IR CIB corresponds to gray- 
body radiation with T ~ 18K, the AGN CIB mainly 
comes from high redshifts. Utilizing the redshift in- 
formation obtained from the cross correlation with 
galaxies, AGN contamination could be localized to 
high z. 

• AGN-late type galaxy cross correlation is much 
weaker than AGN-early type galaxy cross correla- 
tion at arc minute scale (Brown, Boyle, & Webster 
2001) while SFGs strongly correlate with late type 
galaxies. For future CIB survey which can resolve 
this scale, one can distinguish the contribution of 
AGN by cross correlating CIB with early and later 
type galaxies separately. 

3.1.2. The integrated Sachs- Wolfe effect 

The time variation of gravitational potential field intro- 
duces secondary CMB temperature fluctuations (Sachs & 
Wolfe 1967): 
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If the variation of gravitational potential 



(8) 



is caused by 

the deviation from a TO = 1 universe, this effect is called 
the integrated Sachs- Wolfe (ISW) effect. The contribution 
of the ISW effect to CIB intensity is 



hsw = Oisw 



Jfir du V 



9i SW 0.2 nW nU 2 sr" 1 . (9) 



The temperature fluctuation caused by ISW <50isw is 
< 10~ 5 , so fluctuations of the CIB caused by ISW are 
~ 10~ 6 nW m~ 2 sr _1 . We then can safely neglect this 
effect. 

3.1.3. The thermal Sunyaev-Zeldovich effect 

Free electrons scatter off CMB photons by their ther- 
mal motions and introduce secondary CMB temperature 
fluctuations, namely, the thermal Sunyaev-Zeldovich (SZ) 
effect. It contributes an intensity I$z to the CIB (Zel- 
dovich & Sunyaev 1969) 

J Fm exp(x) - 1 V tanh(x/2) J 

~ j/10 3 nW m~ 2 sr _1 . 

Here x = hv/ksTcMB- Since electron thermal energy is 
much higher than CMB photon energy, photons gain en- 
ergy and shift toward infrared frequency after scattering. 
So the SZ contribution to the CIB is much larger than that 
of the ISW effect. 

The y parameter is given by the integral of electron ther- 
mal pressure p e along the line of sight 

f n e k B T e f p e 

JJ= I 2 a Tadx- 



m e c* 



T a T adx- (11) 



The fluctuation of y is <~ 10 6 at degree scale (see, 
e.g., Zhang & Pen (2001); Zhang, Pen, & Wang (2002)). 
Then fluctuations it introduces to the CIB are ~ 
10 -3 nW m~ 2 sr _1 , which arc still ~ 1000 times smaller 
than SCIB fluctuations and can be neglected. 



3.1.4. Brown drawfs 

Brown drawfs associate with galaxies, so the CIB caused 
by brown drawfs has fractional fluctuations of the same 
order as the SCIB. Since the total mass of brown drawfs 
accounts for <C 1% of the total energy in the universe, 
its contribution to CIB intensity is <C 0.01nWm _2 sr _1 
(Karimabadi & Blitz 1984). So its contribution to CIB 
fluctuations is at most ~ 0.1% of SCIB fluctuations and 
can be neglected. 



3.1.5. Intergalactic dust 

Intergalactic dust correlates with galaxies. The CIB 
contributed by intergalactic dust should have similar frac- 
tional fluctuation as the SCIB. But the contribution of 
intergalactic dust to the CIB flux seems to be negligi- 
ble based on SCUBA measurement (see Hauser & Dwek 
(2001) and reference therein for detailed discussion). The 
existence of abundant intergalactic dust would dim far- 
away SNIas and alter the conclusion of accelerating expan- 
sion of the universe based on SNIa Hubble diagram. Since 
WMAP has independently confirmed the acceleration of 
the expansion of the universe without SNIa prior (Spergel 
et al. 2003), there is little room left for intergalactic dust. 
So we will neglect this possible source. 



3.1.6. Other cosmic sources 

Bond, Carr, & Hogan (1986) discussed several other pos- 
sible CIB sources. Primeval galaxies, first stars and decay- 
ing particles all reside at very high redshifts. Gravitational 
lensing caused by low z dark matter (shallower than the 
corresponding galaxy survey) correlate the high z emis- 
sion with low z galaxies. But such correlation is negligible 
due to two reasons. Firstly, the correlation is caused by 
low z dark matter lensing. For example, the cross correla- 
tion with SDSS would be mainly caused by dark matter at 
z z nl /2 = 0.2. So the lensing effect is small. Secondly, 
these high z CIB sources only contribute a small fraction, 
if not negligible, to the total CIB intensity. 



3.1.7. Contaminations of galaxy surveys 

Galactic dust causes extinction, which biases the ob- 
served galaxy distributions. This effect may correlate the 
galactic dust foregrounds of CIB with extragalactic galax- 
ies. But since at least the first order effect of extinction 
is included in the observed galaxy number distribution 
dn/dz, one expects the residual correlation caused by dust 
extinction to be negligible. Stellar contamination in galaxy 
surveys also correlate CIB foregrounds with galaxy distri- 
bution. But such contamination is at 1% level (e.g. SDSS, 
Sheldon et al. (2003)) and is thus negligible, too. 
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3.2. Statistical errors 
The statistical errors of Cir^g measurement are 



A(CJS G ) 



1 + (Cm. + C dust + %¥ t )(C G + 



TR,G 



(2/ + i)A;/ sky 



(12) 

In this expression, we have assumed Gaussian fluctua- 
tions in the CIB, which should at least hold at degree 
or larger scales. In the [0.3,3] THz frequency range, 
the integrated CMB intensity fluctuation is negligible. 
We consider two dominant noise sources: interplanetary 
and galactic dust, and short noise of the signal. We 
adopt Cdust = (10 nWm^sr 1 ) 2 ! -3 . The I dependence 
is adopted from Wright (1998) and the amplitude roughly 
corresponds to the case avoid of the galactic plane (Wright 
1998). The short noise power spectrum is estimated by 
Cshort = 47r/iVsFG- ln this expression, we have omitted 
the flux variation of individual SFG. We further assume 
the total number of SFG to be N SFG = 2 x 10 9 . The 
window function Wi reflects the angular resolution of 
the CIB survey. We approximate it as a Gaussian func- 
tion Wi = exp(-l 2 6%/2) with O = FWHM/(2 V / 21n2). 
While the FIRAS beam FWHM is nominally 7° across it 
has sharper edges than a Gaussian beam so for many 
applications it is better approximated by a Gaussian 
beam with 5° FWHM 3 . C NyG = 47r/ sky /7V G is the Pois- 
son noise of a galaxy survey. We adopt the expected 
number of SDSS galaxies as iV G = 5 x 10 7 and choose 
/sky = 1/4. The bin size Al is chosen to be Al = 0.5Z. 
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Fig. 2. — The statistical errors in the CIB-galaxy correlation 
measurement. The rcdshift bins adopted arc: [0, oo] (solid line), 
[0.0,0.2] (dot line), [0.4,0.6] (short dash line), [1.0,1.2] (long dash 
line) and [1.6, 1.8] (dot-dash line). The top panel corresponds to 
FIRAS+SDSS and the bottom panel corresponds to a toy CIB ex- 
periment with 1° FWHM. The cross correlation measurement is 
mainly limited by the angular resolution of CIB experiment at low 
2 and limited amount of galaxies at high z. 



Cir.G can be measured with <~ 10% accuracy at 10° scales 
and' z < 0.6 for a combination of FIRAS and SDSS (Fig. 
2). In these ranges, the error of cosmic variance dominates. 
To beat down the cosmic variance, higher angular resolu- 
tion CIB experiment is required, which can probe smaller 
angular scales (Fig. 2). Since the median redshift of SDSS 
galaxies is z m <~ 0.5, where most correlation signal should 
come from, one expects the optimal measurement to be 
at z ~ 0.5 (Fig. 2). Since only a small fraction of SDSS 
galaxies lies at > 1, the correlation signal is weak and the 
short noise of galaxy distribution is large. Thus the corre- 
lation measurement becomes noisy. To measure the cross 
correlation at z > 1.5, a deeper galaxy survey is needed. 
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Fig. 3. — The estimated statistical errors of measured from the 
CIB-galaxy correlation. Data points connected with solid lines cor- 
respond to FIRAS and those connected with dash lines correspond 
to a toy CIB experiment with 1° FWHM. Filled circle points corre- 
spond to SDSS with a Galaxy median redshift z m = 0.5 and filled 
triangle points corresponds to a toy galaxy survey with z m = 1.0 
and the same amount of galaxies at low redshift as SDSS. We have 
assumed the 3D correlation of IR emissivity-galaxy to be known. 

Given A 2 - Gl jd{z) can be directly inferred from Cir.g 
(Eq. 6). The accuracy of the measured jd(z) is then de- 
termined by the accuracy of Ctr,,g measurement. We ap- 
proximate the relative error of the inferred jd(z) as the 
minimum ACij^g/Cir.g- The result estimated by this 
somewhat arbitrary simplification should not deviate from 
that of a sophisticated analysis by a factor ~ 2 (Fig. 2), so 
it is sufficient for our purpose. Currently, the extraction 
of jd(z) is mainly limited by the finite depth of the galaxy 
survey (Fig. 3). FIRAS+SDSS could measure the dust 
emissivity jd with ~ 10% accuracy at z < 0.7. For higher 
z, the measurement becomes noisy quickly. The accuracy 
of jd measurement is limited by the angular resolution of 
CIB experiment at low z and the amount of galaxies in 
galaxy surveys at high z. Combining a possible future 
CIB experiment with 1° FWHM and a deeper galaxy sur- 
vey with median redshift z m = 1.0, one may probe jd and 



3 Private communication with Dale J. Fixsen. 
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the star formation history to z ~ 1.5 with several percent 
statistical errors. 

Since our redshift bins are Az ~ 0.2 and one does not 
expect strong evolution of A 2 G over these redshift ranges, 
the above analysis does not require precision measurement 
of galaxy redshifts. Photometric redshift can be measured 
with cr z /(l + z spcc ) <~ 0.04 at z < 6 (Massarotti, Iovino, 
Buzzoni, & Valls-Gabaud 2001), so the error of photo-z 
should not affect our result significantly. 

The extraction of jd is based on known A 2 G , which is 
in principle an observable. Since at sufficiently large scale, 
the cross correlation coefficient of jd and galaxies is unity, 
combining the mean CIB flux and CIB auto correlation 
measurement and galaxy A G measurement, jd and A 2 G 
can be determined simultaneously given the cosmology. 
However, the feasibility of this approach relies on future 
accurate removal of foregrounds, so we postpone its dis- 
cussion in this paper. Furthermore, as we will show in §4, 
the measurement of CIB-ISW cross correlation is able to 
determine the mean bias of jd and thus A 2 G . 

4. CIB-CMB CROSS CORRELATION 

Recent successful detections of the cross correlation be- 
tween WMAP and galaxies (Fosalba & Gaztanaga 2003; 
Fosalba, Gaztanaga, & Castander 2003; Scranton et al. 
2003; Afshordi, Loh, & Strauss 2003) present a power- 
ful way to constrain cosmology, especially dark energy 
(through the ISW effect). But current galaxy surveys are 
either limited by finite sky coverage (For APM, 4300 deg 2 
and for present SDSS, 3400 deg 2 ) or by finite survey depth 
(For 2MASS, (z) < 0.1). This limits current detections to 
< 3<T level. To bypass these limitations, one can cross 
correlate CMB with other cosmic backgrounds, to reduce 
cosmic variance and to amplify cross correlation signal. 
One choice is the X-ray background (XRB). But the mea- 
surement of WMAP and the all sky XRB survey ROSAT 
fails to show any significant correlation (Diego, Silk, & 
Sliwa 2003). On one hand, XRB is dominated by AGNs, 
which mainly concentrate to high redshift. On the other 
hand, ISW mainly comes from low redshift where £l m de- 
viates significantly from 1. At degree scale, which is the 
finniest scale WMAP+ROSAT can probe, SZ contribu- 
tion also mainly comes from low redshift (see, e.g., Zhang 
& Pen (2001)). So one does not expect a strong correla- 
tion between WMAP and ROSAT. Another choice is the 
cosmic radio background (CRB). Since the Milky Way is 
a strong synchrotron radiation source, which has a similar 
spectral index as the CRB synchrotron component, the 
subtraction of CRB foregrounds is difficult. Without a 
robust subtraction of CRB foregrounds, the tight correla- 
tion between galactic foregrounds (For CMB, synchrotron 
and free- free emission, Bennett et al. (2003)) may bias the 
CIB-CMB correlation measurement significantly. On the 
other hand, the origin of CRB is difficult to predict ana- 
lytically. The synchrotron component of CRB associates 
with magnetic field, which is poorly understood. The ther- 
mal component of CRB associates with HII regions gen- 
erated by massive star formation and is hard to predict 
too. So, even if the intrinsic cross correlation is detected, 
it would be hard to interpret. MeV cosmic 7-ray back- 
ground (CGB) is expected to have a tight correlation with 
the large scale structure, but current CGB experiments 



may be too noisy to detect any significant cross correla- 
tion with CMB (Zhang & Beacom 2004). 

SCIB traces star forming galaxies, which have strong 
correlation with gravitational potential (ISW) and gas 
thermal pressure (SZ effect). The modeling of SCIB is 
relatively straightforward, as discussed in previous sec- 
tions. So, the interpretation of CIB-CMB correlation is 
relatively robust. Thus, SCIB-CMB cross correlation is 
a good way to measure the ISW effect and the SZ effect. 
The detection of such cross correlation signal is observa- 
tionally feasible, as implied by the successful detection 
of COBE-FIRAS correlation (Burigana & Popa 1998). 
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Fig. 4.— CIB-CMB cross correlation. The SZ effect is esti- 
mated in the Raylcigh-Jeans regime. The dot lines are the upper 
limit of systematic errors caused by the CIB and the CMB dust 
foreground correlation. This estimation assumes no foreground re- 
duction in CMB maps and assumes perfect foreground correlation. 
With reasonable CMB foreground removal, the systematic error can 
be reduced by a factor of 10. 



4.1. CIB-ISW cross correlation 

The cross correlation power spectrum of ISW-CIB is 
given by 



^C m , ISW I = J ^J±-^ 2 ^^{Lz ul\. (13) 



If assuming a bias model (Sj = bjS), we obtain 



H ^\Jt {a J^ ] - a 2P 



A 2 ;(fc,z) 
3<P y ' ' 2 



(14) 



3„ 2 , d D A 2 
" ' 2 Hobj dt { ^DV ] linearre S imc 
Here, D is the linear density growth factor, calculated by 
the fitting formula of Carroll, Press, & Turner (1992). a is 
the scale factor. H is the present Hubble constant. The 
power of A 2 ^ concentrates on large scale due to the k 2 

term in the denominator, thus the ISW effect dominates 
the CIB-CMB correlation at large angular scale (I < 100) 
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and drops to zero quickly toward small scale (Fig. 4). The 
correlation amplitude reaches ~ 0.5^iK at I ~ 20. One 
bonus of ISW-CIB correlation is that it gives the mean bj 
averaged over redshift. 

4.2. CIB-SZ cross correlation 

The CIB-SZ cross correlation power spectrum in the 
Raylcigh- Jeans regime is given by 

^C IR;SZ / = -2 / kJ^ aT ^A%(^z)d X . (15) 
2tt J A m e c z I n x 

Here, A 2 p is the cross correlation power spectrum (vari- 
ance) of Sj and 5 P = p e /pe — 1. We calculate the gas 
density weighted temperature p e by the continuum field 
model (Zhang & Pen 2001), which agrees with ACDM 
and self similar hydro simulations very well (Zhang, Pen 
& Trac 2004). This model predicts a mean temperature 
decrement 5.0/zK. The modeling of Aj is much more com- 
plicated. For simplicity, we assume S p traces dark matter 
overdensity with a constant bias b p = 5, as predicted by 
Zhang & Pen (2001). The CIB-SZ correlation dominates 
at I > 100. At one degree scale (I ~ 300), the cross corre- 
lation amplitude reaches — 0.3/iK (Fig. 4). 

With a possible future CIB experiment of (sub)degree 
angular resolution, the CIB-SZ cross correlation can be 
measured. The SZ effect is sensitive to various thermal 
heating processes such as SN feedback and quasar feed- 
back. The O and B star formation responsible for SCIB 
is also responsible for SN feedback on the intcrgalactic 
medium, so the measurement of CIB-SZ cross correlation 
helps to constrain the role of SN feedback on the SZ effect. 

4.3. Observational feasibility 

Since we will utilize (almost) full sky data, we expect 
a factor of 3-4 decrease in the sample variance with re- 
spect to current WMAP+galaxy measurement. Further- 
more, Interplanetary dust does not correlate with CMB 
foregrounds. The main obstacle of CIB-CMB cross cor- 
relation measurement comes from the correlations of CIB 
and CMB galactic foregrounds, mainly the galactic dust 
foregrounds. But such correlations concentrate on large 
scales. The power spectrum of CIB galactic dust fore- 
ground scales as C; oc l~ 3 . CMB galactic foregrounds 
have a combined power spectrum C; cx l~ 2 (Bennett et al. 
2003). These foreground correlations degrade the detec- 
tion of the ISW effect, but has only minor effect on the 
detection of the SZ effect. Furthermore, one can avoid the 
galactic plane to minimize foregrounds. 

We can estimate the upper limits of the systematic er- 
rors caused by foreground correlations. For the WMAP W 
band, the only non-negligible foreground is galactic dust 
emission with a flat power spectrum Cil 2 /(2ir) ~ 50yuk 2 
(Fig. 10, Bennett et al. (2003)). This foreground strongly 
correlates with the CIB galactic foreground, which has a 
power spectrum d ~ (10 nWrn^sr" 1 ) 2 /- 3 (§3.2). With- 
out any foreground removal, the systematic error intro- 
duced by foreground correlations prohibits the CIB-CMB 
correlation measurement at I < 20, but it still leaves a 
window at 20 < I < 60 for the CIB-ISW cross correla- 
tion detection. Its effect to the CIB-SZ cross correlation 
measurement at I > 200 is negligible (Fig. 4). 



Current foreground removal method using multi- 
frequency maps can reduce CMB foreground contamina- 
tions in the CMB power spectrum by more than 90%. For 
cleaned WMAP map at W band, foreground only con- 
tributes 0.4% of the CMB power spectrum (Bennett et al. 
2003). Since maps used for foreground removal correlate 
with SCIB, in principle, the foreground removal process 
may introduce false correlation signal with SCIB or lose 
correlation signal. But The ISW contribution to CMB at 
I ~ 10 and the SZ contribution to CMB at I > 300 are 
much larger than 0.4%, thus the CMB signal responsible 
for CMB-CIB correlation do survive after the foreground 
removal process, so we can neglect such complexity in the 
cleaned map. Thus using cleaned CMB map, one expect 
a factor of ~ 10 reduction in the systematic errors. 

FIRAS has an angular resolution ^5° and thus can 
probe I ~ 20. Around such scales, the statistical error of 
CIB-CMB correlation measurement is ~ y/2/(2l + l)Al < 
10% and the systematic error is < 10% using cleaned 
WMAP W band map. So, combining WMAP and FI- 
RAS, one expects ~ 20% accuracy in the CMB-CIB cross 
correlation measurement. 



5. CONCLUSION 

Cosmic infrared background (CIB) traces the large scale 
structure of the universe and contains plenty of informa- 
tion of the star formation history. But overwhelming fore- 
grounds prohibit its precision measurement. By cross cor- 
relating cosmic infrared background with galaxies, one can 
eliminate CIB foregrounds, minimize and localize back- 
ground contaminations and directly obtain redshift infor- 
mation of CIB sources from galaxy photometric redshift. 
Since one does not rely on spectral information to extract 
redshift information of CIB sources, one only needs to esti- 
mate the integrated CIB, which is directly determined by 
star formation rate through the energy conservation and 
is thus fairly model independent. 

The CIB-galaxy cross correlation at degree scale is 
<~ 10% and can be measured with 10% accuracy. We 
estimated that, the cross correlation could enable a di- 
rect and statistically robust measurement of star forma- 
tion rate with 10% to z < 1.5. With our model based on 
the integrated CIB intensity, we predict that the cross cor- 
relation between the CIB and the CMB is about 0.5fiK at 
10° scale (the ISW effect), changes sign at several degree 
scale and reaches ~ — 0.3/xK at one degree scale (the SZ 
effect). We argue that this cross correlation is observa- 
tionally feasible. Such measurement would constrain the 
amount of dark energy by the ISW effect, the amount of 
the thermal energy of the universe by the SZ effect and 
provide further constraint on the clustering property of 
SFG and make our CIB model self consistent. 
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